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TO  A  POISSON  SHADOWING  PROCESS 

By 
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Abstract 

A  star  shaped  curve,  C  ,  in  the  plane  is  subject  to  a  Poisson  shadowing 

process.  According  to  this  process,  disks  of  random  size  appear  at 

random  locations  in  a  region  between  a  source  of  light,  which  is  at  the 

origin,  and  the  curve  C  .  These  disks  cast  shadows  on  C  .  Trials  are 

conducted  along  the  lighted  portion  of  C  .  Each  trial  requires  a  fixed 

length,  Z  ,  of  C  .  The  different  trials  are  independent  and  have  a 

fixed  probability,  p  ,  of  success.  The  number  of  trials  conducted  along 

C  is  a  random  variable,  N  ,  which  depends  on  the  random  length  of  the 

lighted  portion  of  C  .  The  success  probability  is  P  =  1  -  E{q‘ }  »  where 

✓ 

q-l-p  .  Lower  and  upper  bounds  for  P  are  derived.  A  numerical  example 
shows  cases  in  which  these  bounds  could  be  very  close. ^ 
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1.  Introduction 

Consicer  a  star-shaped  curve  in  the  plane,  C  ,  and  a  source  of 
light  at  the  origin,  0  .  If  there  are  no  obstacles  between  the  origin 
and  the  curve  C  ,  the  whole  curve  is  in  the  light  or  visible.  Certain 
experiments  (trials)  can  be  conducted  along  the  lighted  portion  of  the 
curve.  Each  such  trial  requires  a  length  L  of  C  ,  and  the  probabil¬ 
ity  of  its  success  is  p  ,  0<p<l  .  Let  L{C}  be  the  total  length  of 
C  .  If  C  is  completely  visible,  N  -  [L{C}/£]  trials  can  be  con¬ 
ducted.  ([a]  designates  the  integer  part  of  a.)  Assume  that  all  these 
trials  are  independent,  having  the  same  success  probability,  p  , 

(Bernoulli  trials).  Thus,  the  probability  of  at  least  one  success,  when 

N 

C  is  completely  lighted,  is  S  -  1  -  q  ,  where  q  -  1  -  p  .  In  reality, 

C  may  not  be  completely  visible,  due  to  shadows  cast  on  it  by  objects 
(disks,  say),  which  are  randomly  dispersed  in  a  region  between  0  and  C  . 
The  centers  of  the  disks  follow  a  Poisson  process  and  their  diameters  are 
of  random  size.  Thus,  the  total  number  of  trials  that  can  be  performed 
along  the  visible  portion  of  C  is  the  random  variable 
K 

N  »  I  [X./£]  ,  where  K  is  the  random  number  of  connected  subsets 
i-1  1 

(disjoint  segments)  of  C  ,  which  are  visible,  and  is  the  length  of 

the  i-th  such  subset.  The  probability  of  at  least  one  success,  under 

N 

this  random  shadowing,  is  P  -  1  -  Q  ,  where  Q  »  E{q  }  .  The  present 
study  develops  a  method  for  determining  upper  and  lower  bounds  for  Q  . 

This  method  is  based  on  the  methodology  developed  by  Yadin  and  Zacks  [1], 
for  determining  the  moments  of  the  total  visible  portion  of  C  , 
k 

V{C}  ■  l  X.  .We  show  that  a  lower  bound  for  Q  (an  upper  bound  for 
i-1  1 

P  ,  respectively)  is  the  value  of  the  moment  generating  function  (MGF)  of 
V{C}  ,  at  the  point  N  -  £nq ft  .  The  upper  bound  for  Q  (lower  bound  for 

ft  JJft 

P)  can  be  obtained  by  considering  Q  -  E{q  )  ,  where 
k  (\~l\  * 

N  -  l  l— y — I  ,  and  a,  -  max(a,0)  .  In  section  2  we  present  the 

i-1  '  '+  + 

structure  of  the  Poisson  field  of  shadowing  disks  and  review  the  main 


results  of  [1],  which  lead  to  the  MGF  of  V{C}  .  In  Section  3  we  present 
lower  and  upper  bounds  for  Q  ,  which  are  based  on  the  MGF  of  V{C}  and 

*  K 

V  {C}  -  i  (X.--0  . 

i-1  1 

Section  4  presents  a  numerical  example  for  annular  regions  and  a 
standard-uniform  Poisson  field  of  shadowing  disks.  It  is  shown,  in  the 
numerical  example  of  Section  4,  that  the  lower  and  upper  bounds  for  Q 
developed  in  the  present  paper  could  be  very  close  and  effective.  The 
present  paper  is  very  tightly  linked  with  our  previous  paper  [1],  in 
which  the  general  theory  for  the  determination  of  the  moments  of  visi¬ 
bility,  ■  E{Vn{C}}  ,  is  presented.  As  such  it  could  be  considered 
as  an  extension  of  [1]  for  an  important  class  of  applications. 
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2.  The  Poisson  Field  of  Shadowing  Obstacles,  and  the  Moments  of  V{C>  . 

Consider  a  star-shaped  curve,  C  ,  given  by  a  function  r(s)  on  an 
interval  [s',s"]  ,  i.e., 

C  *  {(p.6)  ;  P  *  r(6)  ,  s'  £  0  £  s”)  .  (2.1) 

We  further  assume  that  shadows  on  C  are  cast  by  disks,  which  are 
randomly  distributed  within  a  region,  ,  bounded  by  the  curves 

U  ■  l(p,0)  ;  P  ■  u(9)  ,  s*  <  6  £  s**}  and 

W  *  {(p,0)  ;  p  *  w(0)  ,  s*  £  0  £  s**}  .  Each  disk  is  characterized  by  a 

point  (p,6,y)  ,  where  (p,0)  are  the  polar  coordinates  of  its  center 

and  y  is  its  diameter.  It  is  assumed  that  the  centers  are  uniformly 

distributed  between  U  and  W  and  the  diameters  of  the  disks  are  i.i.d. 
random  variables  having  a  c.d.f.  G(y)  concentrated  on  [a,b]  (the  standard  case). 
Moreover,  u(0)  ,  w(0)  and  b  are  such  that  both  the  origin,  0  ,  is 
uncovered  and  the  curve  C  is  not  intersected  by  any  one  of  the  random 
disks.  For  the  precise  conditions  see  [  1  ].  It  is  further  assumed  that 
the  number  of  disks  whose  centers  fall  within  a  subset  C  of  has  a 

Poisson  distribution,  with  mean  A  ,  when  A  is  the  area  of  C  . 

A  point  P  ■  (r(s),s)  is  said  to  be  visible,  if  the  ray  OP  is  not 
intersected  by  any  shadowing  disk. 

The  measure  of  visibility  V{C}  is  defined  as 

s" 

V(CJ  -  /  I(sK(s)ds  ,  (2.2) 

s' 

where  I(s)»l  if  the  point  (r(s),s)  is  in  the  light  (visible),  and 
l(s)«0  otherwise.  £(s)ds  -  (r2(s)  +  (r ' (s))2]1/2ds  is  the  infinitis- 
oal  length  of  C  at  (r(s),s)  .  The  moments  of  V{C)  where  expressed 
ia  [  1]  in  terms  of  the  K-functions,  K_(s,t)  and  K+(s,t)  .  yK_(s,t) 
and  yK+(s,t)  are,  respectively,  the  expected  number  of  disks  in  , 
whose  centers  have  orientation  coordinates  in  [s-t,s]  ([s,s+t],  resp.), 
and  which  do  not  Intersect  the  ray  with  orientation  s  .  It  is  shown  in 
[  1]  that  these  functions  are  given  by 
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K  (s,t) 


and 


K+(s,t) 


8  w(6) 

/  /  G(y(p ,8-6) )pdpd8 
s-t  u(9) 


s+t  w(0) 

/  /  G(y (p,6-s))pdpd9 

s  u(0) 


where 


y(p,9-s) :  y(p,s-9) 


2psln js-9 | 


,  if  I  s-6  |  <tt/2 


,  if  ls-9|27r/2 


is  the  maximal  diameter  of  a  disk  centered  at  (p,9)  ,  which  does  not 
intersect  the  ray  with  orientation  s  . 

It  is  shown  in  [  1 J  that  the  n-th  moment  of  V(C)  is 


(2.3) 


(2.4) 


P  *  n!  /.  .  .  / 
n 


p(s.  ,...,s  )  II  Z(s  )ds .  i 

s'*s.<;...ss  s8"  1  n  i-l  1  1 

1  n 


(2.5) 


where  p(s^, . . . ,8^)  is  the  probability  that  n  points  on  C  ,  with 
orientation  coordinates  si»,,,»sn  are  simultaneously  visible.  It  Is 
further  shown  that 

p(s1,...,8n)-  exp^vtCj^ exp  ■jyK_  (S^S^S*)  + 


(2.6) 


VK. 

in  which 


+“  "f*  (w  ^T"1)  +  S-C=1+r  j 


I 
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v{G1> 


s** 

-B  / 

2  * 
s* 


^2(s)-u2(s)jds 


Furthermore,  let 


^(s)  ■  exp{yK_(8,s-sr)}  , 

and  define  recursively,  for  j  s  1 
s 

¥j(s)  -  I  ^  ^CyJF^Cy)  exp{y  K_(y,^jZ)  +  UK+Cy,^^)}dy 


yn  -  n! 


^(s)<Pn-1(s)  exp{pK+(s,s**-r.)}ds 


Then 
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3.  The  MGF  of  V{C>  and  The  Bounds  for  Q 

In  Section  1  we  introduced  the  random  variables  K.X^^j  . . .  ,X^  , 
which  are  the  number  of  lighted  (visible)  disjoint  segments  of  C  , 
and  their  length.  Accordingly, 

V{C}  -  I  X  .  We  also  defined  the  random  variable  N  -  I 
i-l  1  i-1 

Thus  N  s  V{C}/£  ,  with  probability  one.  It  follows,  for  every 

q,  0<q<l  ,  that. 

Q  -  E{qN}  ^  E^V{C}/^«  ^(lnqM)  ,  (3.1) 

where  M^(u)  is  the  MGF  of  V(C}  .  Thus,  My  (lng/Z)  is  a  lower  bound 
for  Q  .  Notice  that,  since  V(C)  5  L(C)  <  ®  ,  all  the  moments  of 
V(C)  are  bounded  by  powers  of  1(C)  .  Hence,  the  MGF  of  V(C)  can  be 
expressed  as 

“  i 

M(,(u)  »  z  TT  y.  ,  -“<u<°°  .  (3.2) 

i-0  11  1 


For  the  derivation  of  an  upper  bound  for  Q  ,  we  consider  the  random 
k  X  -l 

variable  N*  -  I  — j —  ,.  Since  N*  <  N  with  probability  one, 

i-l 


In  order  to  obtain  Q*  we  define  a  new  visibility  measure 

.  fe 

V*(C)  -  I  (X,-£).  (3.4) 

i-l  x 


The  moments  of  V*(C)  can  be  obtained  by  the  formulae  presented  in 
Section  2,  in  which  the  K-f unctions  in  (2.6)  -  (2.10)  are  modified  in 
the  following  manner.  Replace  K_(s,t)  and  K+(s,t)  by  K_(s-t1(s)  , 
(t-T^C*))^)  and  K+(s+Tj(s)  ,  (t-t^a))^  ,  respectively,  where  t1(s) 
should  be  determined  so  that 


/  £(s)ds  -  II 2 

•-Ti(s) 


(3.5) 
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and,  similarly,  t2(s)  should  satisfy  the  equation 


s+t2(s) 

/  £(s)ds  -  ^/2  .  (3.6) 

s 

By  definition,  K+(s,0)  -  0  for  all  s  . 

ff 

More  specifically,  let  (n=l,2,...)  be  the  n-th  moment  of  V*(C)  , 
which  is  given  by 


* 

y  *  n If  ...  f  p*(s1,...,s  ) 
n  gin 

where  S-vs'  £s-S...£s  s 

l  1  * 

The  function  p*(s^, . . ..s^)  is  the  probability  that  the  union  of  n 
segments  of  C  ,  each  one  of  length  t  ,  centered  around  the  points 
(r(s^),Sj)  ,  i-l,...,n  ,  is  completely  visible.  v.n  other  words,  define 
the  indicator  function  I*(s)  ,  which  is  equal  to  1  if  the  segments 
of  C  of  length  t  ,  centered  at  (r(s),s)  ,  is  completely  visible  and 
is  equal  to  0  otherwise.  Accordingly, 


n  £(s.)ds,  ,  (3.7) 

1-1 


» 


V*(C)  -  /  I^(s)Z(s)ds 


(3.8) 


As  explained  above,  p*(s1,...,s  )  -  E{  n  I  (s.)} 

X  n  i=l  1 


Following  the  theory 


developed  in  [  i ]  ,  p*(s^, . . . ,sn)  is  given,  as  in  (2.6),  by 
p*(s1,...,sn)  -  exp^-v{C1^  exp|yK_^1-T1(s1)  , 


(3.9) 


Finally,  if  MyA(u)  denotes  the  MGF  of  V*(C)  ,  then  Q*  ■  M^Cln  q/£)  , 
which  is  the  upper  bound  for  Q  .  In  the  following  section  we  illustrate 
these  bounds  in  a  special  case. 
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4.  Lower  And  Upper  Bounds  For 


in  A  Special  Case 


In  the  present  section  ve  exhibit  the  method  of  determining  lower 
and  upper  bounds  for  the  failure  probability  Q  .  in  the  following  special 
case.  The  curve  C  is  an  arc  on  a  circle  of  radius  r  ,  centered  at  the 
origin,  limited  by  rays  having  orientations  s'  and  s"  , 

TT  IT 

-  ~2  <s*  £  s"  s-2*  centers  the  disks  are  distributed  between 

U  and  W  ,  where  U  and  W  are  circles  centered  at  the  origin,  with 
radii  0  <  u  <  w  <  r  .  Moreover,  we  assume  that  the  centers  of  the 
disks  are  uniformly  distributed  within  this  annular  region,  and  their 
random  diameters,  Y  ,  are  uniformly  distributed  between  [a,b]  ind 
pendently  of  their  centers,  where  0<|su<wir-|.  This  spe« 
case  was  previously  studied  in  [  1  ] .  We  have  shown  that  in  the  pres 

A  A 

case,  K_(s,t)  -  K+(s,t)  =  K(t)  .  Explicit  formulae  for  this  functii 
can  be  found  in  [  1  ] .  Notice  that  in  the  present  case  of  C  being  a 
circular  arc,  t. (s)  -  t0(s)  -  t/2x  for  all  s  .  Accordingly,  in  the 
determination  of  Q*  we  replace  K(t)  by  K  (t  -  ^)+  .  The  compu¬ 

tation  of  the  moments  u  and  y*  follows  the  procedure  described  in 

n  n 

[1]. 

Let  {r)n  ,  nil}  be  the  normalized  moments  of  V{C}  ,  i.e., 

nn  ■  E{Vn(C)/rn(s' '-s')n}  .  The  sequence  (nn  ;  nil}  is  decreasing  and, 

as  shown  in  [  1  ]  ,  lim  n  »  P.  ,  which  is  the  probability  of  complete 

n-*°  n 

»  „n 

visibility  of  C  .  Furthermore,  My(lnq/£)  -  Z  nQ  ,  where 

0  ■  (lnq)  r(s"-s')/£  .  One  can  approximate  this  MGF,  to  any  degree  of 
accuracy,  in  the  following  manner.  Given  an  arbitrary  g  ,  g>0  ,  let 
m  be  a  positive  integer  such  that 


(4.1) 


1 

I 
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Let 


Then,  S 


<W-1S 


1,  + 

i  m 


ra-i  0i 

ex  -  Z  jr 
1=0 


(4.2) 


Indeed, 


■  ii>  (•■-•.)  Hi, 
*  (■•  -  '■)  (• 


9 1 

|e  1  -  Z 


(4.3) 


can  be  approximated  in  the  same  manner.  In  Table  1  we  provide 
numerical  values  of  the  lower  and  upper  bounds  for  Q  ,  corresponding  to 
the  following  parameters  of  an  annular  region:  r=1.0  ,  w=.6  ,  u=.4  . 

The  parameters  of  the  distribution  of  Y  are  a=.l  and  b=.5  .  These 
bounds  are  given  for  two  values  of  p  ,  two  values  of  £  *  two  values 

of  A  =  s''  -  s'  ,  and  q  =  .8  . 

Table  1.  Lower  and  Upper  Bounds  for  Q  ,  Circular  Arc,  C  ,  And  Annular 
Region  of  Disk  Centers 


1 

A  =  120° 

> 

II 

o 

r. 

P=1 

P=5 

y=l 

P=5 

.£=.2 

.126861 

.117760 

.258360 

.210320 

.353190 

.341425 

.502616 

.454504 

£=.4 

.356616 

.337353 

.521052 

.439304 

.596354 

.580451 

.721755 

.661923 

This  table  shows  that  in  the  present  case  the  method  developed  here  is 
very  satisfactory. 
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20.  ABSTRACT  (continued) 

i 

on  the  random  length  of  the  lighted  portion  of  C.  The  success  probability  is 
N 

P*1  -  E{q  }  ,  where  q*l-p.  Lower  and  upper  bounds  for  P  are  derived.  A 
numerical  example  shows  cases  in  which  these  bounds  could  be  very  close. 
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